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ABSTRACT 

We use the Sedov-Taylor self-similar solution to model the radio emission from the 7-ray bursts (GRBs) 980703 
and 970508, when the blastwave has decelerated to non-relativistic velocities. This approach allows us to infer the 
energy independent of jet collimation. We find that for GRB 980703 the kinetic energy at the time of the transition 
to non-relativistic evolution, r NR » 40 d, is E ST w (1 -6) x 10 51 erg. For GRB 970508 we find £ ST ~ 3 x 10 51 erg 
at ?nr f» 100 d, nearly an order of magnitude higher than the energy derived in Frail, Waxman & Kulkarni (2000). 
This is due primarily to revised cosmological parameters and partly to the maximum likelihood fit we use here. 
Taking into account radiative losses prior to t^R, the inferred energies agree well with those derived from the early, 
relativistic evolution of the afterglow. Thus, the analysis presented here provides a robust, geometry-independent 
confirmation that the energy scale of cosmological GRBs is about 5 x 10 51 erg, and additionally shows that the 
central engine in these two bursts did not produce a significant amount of energy in mildly relativistic ejecta at late 
time. Furthermore, a comparison to the prompt energy release reveals a wide dispersion in the 7-ray efficiency, 
strengthening our growing understanding that E 1 is a not a reliable proxy for the total energy. 

Subject headings: gamma-rays :bursts — radiation mechanisms:nonthermal — shock waves 



1. INTRODUCTION 

The two fundamental quantities in explosive phenomena are 
the kinetic energy, Eg, and the mass of the explosion ejecta, 
M e j, or equivalently the expansion velocity, [3 = v/c, or Lorentz 
factor, T = (1 -/3 2 )" 1 ' 2 . Together, these gross parameters deter- 
mine the appearance and evolution of the resulting explosion. 
Gamma-ray bursts (GRBs) are distinguished by a highly rela- 
tivistic initial velocity, To J> 100, as inferred from their non- 
thermal prompt emission (Goodman 1986; Paczynski 1986). 
For the range of 7-ray isotropic-equivalent energies observed 
in GRBs, E 7iiso - 10 51 - 10 54 erg (Bloom, Frail & Sari 2001), 
this indicates M e j ~ 10~ 5 - 10~ 3 M Q , compared to several M Q 
in supernovae (SNe). 

The true energy release of GRBs depends sensitively on the 
geometry of the explosion. For a collimated outflow ("jet") with 
a half-opening angle 9j, it is E = fbE lso , where fi, = [1 -cos(6* ; )] 
is the beaming fraction; the true ejecta mass is also a factor of 
fi, lower. Over the past several years there has been growing ev- 
idence for such collimated outflows coming mainly from achro- 
matic breaks in the afterglow light curves (e.g., Kulkarni et al. 
1999; Stanekef al. 1999). The epoch at which the break occurs, 
tj, corresponds to the time at which the ejecta bulk Lorentz fac- 
tor decreases below f?y' (Rhoads 1999; Sari, Piran & Halpern 
1999). 

In this context, several studies have shown that the beaming- 
corrected energies of most GRBs, in both the prompt 7-rays 
and afterglow phase, are of the order of 10 51 erg (Frail et al. 
2001; Panaitescu & Kumar 2002; Berger, Kulkarni & Frail 
2003; Bloom, Frail & Kulkarni 2003; Yost et al. 2003). The 
various analyses are sensitive to the energy contained in ejecta 
with different velocities, V ;> 100 in the 7-rays, T ;> 10 in the 



early X-rays, and T J> few in the broad-band afterglow. How- 
ever, none are capable of tracing the existence and energy of 
non-relativistic ejecta. 

Frail, Waxman & Kulkarni (2000) overcame this problem in 
the case of GRB 970508 by modeling the afterglow radio emis- 
sion in the non-relativistic phase, thus inferring £fK5x 10 50 
erg. This analysis has two significant advantages. First and 
foremost it is independent of jet collimation since the blastwave 
approaches spherical symmetry on the same timescale that it 
becomes non-relativistic (Livio & Waxman 2000). Second, this 
analysis relies on the simple and well-understood Sedov-Taylor 
dynamics of spherical blastwaves, as opposed to the hydrody- 
namics of spreading relativistic jets. In addition, the peak of 
the synchrotron spectrum on the relevant timescale lies in the 
radio band where the afterglow is observable for several hun- 
dred days. 

Two recent developments make similar analyses crucial. 
We now recognize that some GRBs are dominated by mildly 
relativistic ejecta (Berger et al. 2003b). For example, for 
GRB 030329 the kinetic energy inferred from the afterglow 
emission, E K (T ~ few) « 5 x 10 50 erg (Berger et al. 2003b), 
was an order of magnitude higher than the 7-ray energy release 
(Price et al. 2003). Similarly, for GRB 980425 E 1 w 8 x 10 47 
erg (Galama et al. 1998; Pian et al. 2000) was about 1% 
of the relativistic kinetic energy of the associated SN 1998bw, 
E K « 10 50 erg (Kulkarni et al. 1998; Li & Chevalier 1999). 
This begs the question, is there even more energy emerging 
from the engine, either at the time of the burst or later on, at 
non-relativistic velocities? 

Second, there is a growing interest in "unification models" 
for GRBs, X-ray flashes (XRFs) and core-collapse SNe of type 
Ib/c, relying primarily on energetics arguments. For example, 
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namical evolution of the blastwave is described by the Sedov- 
Taylor self-similar solution (Sedov 1946; von Neumann 1947; 
Taylor 1950). In this case the radius of the shock is given 
by r oc (E ST t 2 /A) 1 /^, with p=Ar s . Thus, in the ISM case 
r oc {Esvt 2 /nnip) 1 ! 5 , while in the Wind case r oc (Esjt 2 / 'A) 1 / 3 . 
The constant of proportionality, £(7), depends on the adiabatic 
index of the gas, 7, and is equal to 1 .05 in the ISM case and 
0.65 in the Wind case for 7 = 13/9. The latter is appropriate 
for pressure equilibrium between relativistic electrons and non- 
relativistic protons 3 (Frail, Waxman & Kulkarni 2000). The cir- 
cumburst material shocked by the blastwave is confined down- 
stream to a thin shell of width r/77, with 77 10. 

To calculate the synchrotron emission emerging from this 
shock-heated material we make the usual assumptions. First, 
the relativistic electrons are assumed to obey a power-law dis- 
tribution, N(-f) oc 7" p for 7 > 7 m . Second, the energy densities 
in the magnetic field and electrons are assumed to be a non- 
varying fraction (e^ and e e , respectively) of the shock energy 
density. Coupled with the synchrotron emissivity and taking 
into account self-absorption, the flux received by an observer 
at frequency v and time t is given by (e.g., Frail, Waxman & 
Kulkarni 2000): 

F v =F a {t/t^[(\+z)vf l2 (\-e- T )f,{v/v m )tf{v/v m ), (1) 



the optical depth is given by: 

T„ = T (f /fo) QT [(1 +Z)ls]- (P+4)/2 f2(v/v,n) 1 



and the function 



Lamb, Donaghy & Graziani (2004) argue that GRBs and XRFs 
share an energy scale of ~ 10 49 erg, and that all type Ib/c SNe 
give rise to GRBs or XRFs. Both conclusions result from sig- 
nificantly smaller values of Oj compared to those inferred in the 
past, such that the energy scale, oc 8 2 , is lower by a factor of 
^100 and the true GRB rate, oc 6~ 2 , matches locally the type 
Ib/c SN rate. Given the important ramifications of the GRB 
energy scale for progenitor scenarios we would like to inde- 
pendently address the question: Is the energy scale of cosmic 
explosions 10 49 erg, implicating all type Ib/c SNe in the pro- 
duction of GRBs, or does it cluster on ~ 10 51 erg/ 

The answer will also provide an independent confirmation of 
the jet paradigm by comparison to the isotropic-equivalent ener- 
gies. This is crucial since other explanations for the light curve 
breaks have been suggested, including changes in the density of 
the circumburst medium, a transition to a non-relativistic evo- 
lution on the timescale of a few days (due to a high circumburst 
density), and changes in the energy spectrum of the radiating 
electrons (Dai & Lu 2001; Panaitescu 2001; Wei & Lu 2002). 

Here we address the possibility of significant contribution 
from non-relativistic ejecta and robustly determine the energy 
scale of GRBs independent of geometrical assumptions, us- 
ing Very Large Array 1 radio observations of the afterglows of 
GRBs 970508 and 980703 in the non-relativistic phase. We 
generally follow the treatment of Frail, Waxman & Kulkarni 
(2000), but unlike these authors we carry out a full least-squares 
fit to the data. 

2. THE NON-RELATIVISTIC BLASTWAVE AND FIREBALL 
CALORIMETRY 

The dynamical evolution of an ultra-relativistic blastwave ex- 
panding in a uniform medium (hereafter, ISM) is described in 
terms of its Lorentz factor, T = (\lE KO j'i> r Knm p c 2 r > ~) x l 2 , where 
r is the radius of the blastwave and n is the number density 
of the circumburst medium (Blandford & McKee 1976). This, 
along with the relation for the observer time, which for the line 
of sight to the center of the blastwave is t « r/8T 2 c (e.g., Sari 
1997), determines the evolution of the radius and Lorentz fac- 
tor. For a spherical blastwave the expansion will eventually 
become non-relativistic on a timescale 2 , f NR « 65(£ , ; so ,52/no) 1 ' 3 
d, determined by the condition that the mass swept up by the 
blastwave, M sw w Ei S0 /c 2 , where M sw . 

An initially collimated outflow becomes non-relativistic 

at f NR w 40(£ iso .52/no) 1/4 f],d 4 d ( Livio & Waxman 2000). 
Moreover, as the jet expands sideways (at t ;> tj) the out- 
flow approaches spherical symmetry on a timescale, t s « 

150(£'i so ,52/«o) 1/ ' 4 fj d, similar to f NR . Thus, regardless of the 
initial geometry of the outflow the non-relativistic expansion is 
well-approximated as a spherical outflow. We note that this dis- 
cussion can be generalized to a range of radial density profiles. 
Here, in addition to the ISM model, we focus on a density pro- 
file, p =Ar~ 2 (hereafter, Wind), appropriate for mass loss with 
a constant rate, M w , and speed, v w (Chevalier & Li 2000). 
Following the transition to non-relativistic expansion, the dy- 

1 The VLA is operated by the National Radio Astronomy Observatory, a facility of the National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 

2 Here and throughout the paper we use the notation q = Wq x - 

3 The relative pressure between the protons and relativistic electrons depends on the fraction of energy in relativistic electrons, e e . If this fraction is low, the pressure 
may be dominated by the non-relativistic protons in which case 7 = 5/3. As we show below, e e for both GRB 980703 and GRB 970508 is in the range ~ 0.1 to 0.5 
and thus 7 = 13/9 is applicable. 



/,(*)= / F(y)yW 2 dy. 



(2) 



(3) 



Here, v m = vo(t /to) am /(l +z) is the synchrotron peak frequency 
corresponding to electrons with 7 = 7,,,, F(y) is given in e.g., 
Rybicki & Lightman (1979), and the temporal indices a F , a T 
and a m are determined by the density profile of the circum- 
burst medium. In the 75Mcase a F = 11/10, a T = l-3p/2, and 
ot m = _ 3, while in the Wind case a F = 11/6, a T = — 1— 7p/6, and 
oi m = -7/3 (Waxman 2004a). Equations 1-3 include the appro- 
priate redshift transformations to the rest-frame of the burst. 

Based on the temporal scalings the synchrotron flux in the 
optically-thin regime (v ^> v m ,v a ) evolves as F v oc f( 21_15 rt/io 
(ISM) or F v oc t (5 ~ lp) l 6 {Wind); here the synchrotron self- 
absorption frequency, v a , is defined by the condition T v (v a ) = l. 
Thus, for v^>v m ,v a the transition to non-relativistic expansion 
is manifested as a steepening of the light curves at f NR if the 
outflow is spherical (Sari, Piran & Narayan 1998; Chevalier & 
Li 2000), or a flattening if the outflow was initially collimated 
(Sari, Piran & Halpern 1999). Below, we use this behavior to 
estimate f NR for GRBs 980703 and 970508/ 

In §3 and §4 we use the temporal decay indices and Equa- 
tions 1-3 to carry out a least-squares fit to the data at t > /nr 
with the free parameters Fq, tq, vq and p. These parameters 
are in turn used to calculate the physical parameters of interest, 
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namely r, n e , j m and B; n e m (r)/3)n is the shocked electron 
density (Frail, Waxman & Kulkarni 2000). Since only three 
spectral parameters are available, this leaves the radius uncon- 
strained and thus, 

B = 1 1 J(p + 2)-% 2 _ 52 (r 17 /d 2i ) 4 G, (4) 

7 m = 6J(p + 2)F Q - 52 ^(rn/d2s)~ 2 , (5) 

n e = 3.6 x l0 1 % m Fl_ 52 ^ p)/ \ . 32 r^(r 11 /d 2& r 6 cm" 3 , (6) 

c„ = (1.67 x 10V(5.4 x \0 2 ) (1 ~ p)/2 (p + 2) 2 /(p- 1). (7) 

In the Wind model, the density is appropriate at rsT = r(f NR ), 
i.e., p(r) = nm p (r/r m T 2 . 

To determine the radius of the blastwave a further constraint 
is needed. We note that the energy contained in the electrons 
and magnetic field cannot exceed the thermal energy of the 
Sedov-Taylor blastwave, which accounts for about half of the 
total energy (Frail, Waxman & Kulkarni 2000). The energy 
in the electrons is given by E e = [(p- l)/(p-2)]n e ^ m m e c 2 V, 
while the energy in the magnetic field is E B = B 2 V /Sir; here 
V = 4nr 3 /rj is the volume of the synchrotron emitting shell. 
Thus, using Equations 4-7 and the condition E e + E B < Est/2 
we can constrain the range of allowed values of r. In the ISM 
model Est = nm p (r/\.Q5) 5 [t^/(\ + z)]~ 2 , while in the Wind 
model E ST = A(r/0.65) 3 [f NR /(l + z)T 2 . 

With a constraint on the radius we can also ensure self- 
consistency by calculating the velocity of the blastwave when 
it enters the Sedov-Taylor phase, vst = 2r(l+z)/5fN R (ISM) 
or vst = 2r(l +z)/3/nr (Wind). We expect that roughly v ~ c. 
Finally, the isotropic-equivalent mass of the ejecta is given by 
M e j = \imm p r\ T (ISM) or \kAyst (Wind). The actual ejecta 
mass is reduced by a factor relative to this value. 

3. GRB 980703 

In Figure 1 we plot the radio light curves of GRB 980703. 
The data are taken from Berger, Kulkarni & Frail (2001) and 
Frail et al. (2003). Two gross changes in the light curves evo- 
lution are evident: a flattening at t w 40 d at 4.9 and 8.5 GHz 
and a transition to a constant flux density at late time. The latter 
is due to radio emission from the host galaxy of GRB 980703 
with flux densities at 1.4, 4.9 and 8.5 GHz of 65, 50 and 40 
/iJy, respectively (Berger, Kulkarni & Frail 2001). The flatten- 
ing at t w 40 d marks the transition to non-relativistic evolution 
following a period of sideways expansion of the initially colli- 
mated outflow (Figure 1). A similar value, ?nr ~ 30-50 d has 
been inferred by Frail et al. (2003) from tracking the evolution 
of the blastwave Lorentz factor in the relativistic phase. We 
therefore use here ?nr = 40 d. 

We follow the method outlined in §2 using both the ISM and 
Wind cases. The results of both fits, shown in Figure 1, are 
overall indistinguishable. In what follows we quote the results 
of the ISM model. The best-fit parameters (x 2 min = 123 for 45 
degrees of freedom) are: Fo-52 ~ 2.7, 70,32 ~ 80, 1/0,9 ~ 4.6 and 
p « 2.8. The relatively large values of x^ in is due primarily to 
fluctuations induced by interstellar scintillation, particularly at 
4.9 GHz. 

Using d 2S = dun/Q. +z) 1/2 = 1.4 (z = 0.966, H = 71 km s" 1 

4 We do not consider the Wind case since in this model the observed decay rates 
can be avoided by assuming a break in the electron energy distribution at -y b > 7, 
either 7/, or q. 



Mpc , Cl m = 0.27 and Cl\ = 0.73), and Equations 4-7 we find 
B«1.8x 10~ 2 r 4 l7 G, 7 m w 300r7 7 2 , and n e « 4.9 x 10 3 r7j cm" 3 . 
From these parameters we calculate E c k 3.4 x 10 51 r^ erg, 
£j« 1.7x 10 46 rji erg, and E S t ~ 6.2 x lO 51 ^. These results 
are plotted in Figures 2 and 3. 

The range of blastwave radii allowed by the constraint E e + 
E B <; Est/2 is r i7 w 1.05-2.5, resulting in a range of val- 
ues for the Sedov-Taylor energy, Est ~ (1-6) x 10 51 erg. 
Given the strong dependence on radius, the ratio of energy in 
the electrons to the energy in the magnetic field ranges from 
e e /tB ~ 0.03-9 x 10 4 , while the specific values range from 
e e w 0.01-0.45 and e 8 «5x 10" 6 -0.4. The circumburst den- 
sity is in the range n « 8-3.5 x 10 3 cm -3 , while the blastwave 
velocity is (3st ~ 0.8- 1.9. Finally, the isotropic-equivalent 
mass of the ejecta ranges from (1 -40) x 10" 4 M Q . 

A comparison to the values derived by Frail et al. (2003) 
using modeling of the afterglow emission in the relativistic 
phase is useful. These authors find n w 30 cm" 3 , t e w 0.27 and 
fj-2x 10" 3 . Using the same density in our model (Figure 3), 
as required by the ISM density profile, gives a radius r [7 rs 1.75 
and hence t e w 0.06 and e B ~ 4 x 10~ 3 , in rough agreement; the 
energy is Est ~ 2 x 10 51 erg. 

If we assume alternatively that the energy in relativistic elec- 
trons and the magnetic field are in equipartition, we find r i7 w 
2.05. In this case, Est w 1.5 x 10 51 erg, n w 10 cm" 3 , B w 0.3 
G, and e e = e B = 0.03. 

4. GRB 970508 

The non-relativistic evolution of GRB 970508 was studied 
by Frail, Waxman & Kulkarni (2000). These authors provide 
a rough model for the radio emission beyond r NR w 100 d and 
argue that the constraint E e + E B <J Est/ 2 requires the electron 
and magnetic field energy to be in equipartition, e e = t B K. 0.25, 
with Est ~ 4.4 x 10 50 erg. Here we perform a full least-squares 
fit, using f NR = 100 d, and find somewhat different results. We 
use f NR w 100 d, noting that for GRB 970508 the outflow ap- 
pears to be weakly-collimated (Yost et al. 2003), and hence the 
transition is manifested as a mild steepening of the light curves 
(see §2). 

The best-fit parameters in the ISM model 4 (Xmin =164 for 58 
degrees of freedom) are: Fq- 52 w 38, T0.32 w 3.1 x 10" 3 , ^0,9 ~ 3 
and p Ri 2.17. The large value of x^ in is due primarily to inter- 
stellar scintillation. 

In comparison, Frail, Waxman & Kulkarni (2000) use 
•Pb-52 ~ 41, to.32 ~ 5.3 x 10" 3 , ~ 9.5, and they set p = 2.2; 
a solution with 1/0.9 ~ 4.2 is also advocated but it is not used to 
derive the physical parameters of the blastwave. The formal \ 2 
values for these solutions are 225 and 254, respectively, some- 
what worse than the solution found here. 

As a result, we find that solutions away from equipartition are 
allowed. Adopting the cosmological parameters used by Frail, 
Waxman & Kulkarni (2000), H = 70 km s" 1 Mpc" 1 , Q m = 1 
and Oa = 0, we find Est ~ (6— 11) x 10 50 erg, a factor of about 
20- 100% higher than the values inferred by these authors. 

Using the currently-favored cosmology (§3), we find instead 
that the distance to the burst is higher by about 30%, d 2 $ = 1 .21 
compared to 0.94 (Frail, Waxman & Kulkarni 2000). The 

at 4.9 and 8.5 GHz, F v oc r l 2 , require p ss 1.7 and hence an infinite energy. This 
with a power law index q > 2, but we do not have the data required to constrain 
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change in distance has a significant effect on the derived pa- 
rameters since E e oc eft, E B oc d~ & and Est °c n oc d 6 . Thus, we 
find that the constraint on E e + Es indicates r\i w 3.7-5.9 and 
therefore, fi « 0.04-0.25 G, 7 w 65 - 165 and n « 0.4- 10 
cm" 3 . The Sedov-Taylor energy is £Vr ~ (1.5-3.8) x 10 51 erg, 
while e e « 0.07-0.5 and e B w 0.001-0.45 (Figures 5 and 6). 
Assuming equipartition, we find n 7 = 5.3, E ST = 1.8 x 10 51 erg, 
and e e = efc = 0.11. The derived energy is about a factor of four 
higher than the previous estimate (Frail, Waxman & Kulkarni 
2000). 

A comparison of our best-fit model with the flux of the af- 
terglow in the optical /?-band at t = 110 d, F VjR w 0.3 /zJy 
(Garcia et al. 1998), indicates a break in the spectrum. If 
we interpret this break as due to the synchrotron cooling fre- 
quency, above which the spectrum is given by F v oc v~ p l 2 , we 
find v c w 6 x 10 13 Hz. Since v c = 1.9 x 10 10 fi- 3 (f/110d)- 2 Hz 
we infer B w 0.073 G and hence r xl = 4.3, E ST = 2.8 x 10 5 \ 
e e = 0.25 and t 8 = 8x 10" 3 . These values are in rough agree- 
ment with those inferred from modeling of the relativistic phase 
(Panaitescu & Kumar 2002; Yost et al. 2003), although our 
value of es is somewhat lower. 

5. RADIATIVE CORRECTIONS 

The energies derived in §3 and §4 are in fact lower limits on 
the initial kinetic energy of the blastwave due to synchrotron 
radiative losses. These play a role primarily in the fast-cooling 
regime (v c <C v m ), which dominates in the early stages of the 
afterglow evolution (e.g., Sari, Piran & Narayan 1998). 

Yost et al. (2003) estimate the time at which fast-cooling 
ends, t cm « 0.1 and 1.4 days after the burst for GRB 970508 
and GRB 980703, respectively. Using these values, and our 
best estimate of e e w 0.06 (980703) and e e w 0.25 (970508), 
we calculate the radiative corrections, E oc t m , going back from 
/nr to about 90 s after the burst. Here m « — 17e/12, with 
e = e e /(l + l.Q5e e ) for t < f cm and it is quenched by a fac- 
tor (y m jv c f p ~' l) l' 1 < 1 at later times. Thus, at low values of 
e e the radiative losses are negligible. The cutoff at 90 s cor- 
responds to the approximate deceleration time of the ejecta, 
f dec « 90(£ 5 2/nor^) 1 / 3 s. 

We find that approximately 50% and 90% of the energy was 
radiated away before f NR for GRBs 980703 and 970508, re- 
spectively. Thus, the initial kinetic energies are estimated to 
be 4 x 10 51 erg and 3 x 10 52 erg, respectively. The correc- 
tions from f NR back to t cm , 10% for GRB 980703 and 70% for 
GRB 970508, indicate E K « 2 x 10 51 and 9 x 10 51 erg, respec- 
tively. Both estimates of the energy are in excellent agreement 
with those inferred from the relativistic evolution of the fire- 
ball at f cm (Yost et al. 2003), E K « 3 x 10 51 erg (980703) and 
E K « 1.2 x 10 52 erg (970508). 

6. DISCUSSION AND CONCLUSIONS 

Analysis of the synchrotron emission from a GRB blastwave 
in the non-relativistic phase has the advantage that it is inde- 
pendent of geometry and is described by the well-understood 
Sedov-Taylor self-similar solution. Using this approach to 
model the late-time radio emission from GRBs 980703 (/ > 40 
d) and 970508 (f > 100 d) we infer kinetic energies in the range 
(1-6) x 10 51 erg and (1.5-4) x 10 51 erg, respectively. In- 
cluding the effect of radiative losses starting at fd ec ~ 90 s, we 
find that the initial kinetic energies were about 4 x 10 51 erg and 



3 x 10 erg, respectively. 

The inferred kinetic energies confirm, independent of any as- 
sumptions about the existence or opening angles of jets, that 
the energy scale of GRBs is <~ 5 x 10 51 erg. We therefore un- 
ambiguously rule out the recent claim of Lamb, Donaghy & 
Graziani (2004) that the energy scale of GRBs is of the order 
of 10 49 erg. Since the claimed low energies were based on the 
apparent correlation between E-y^o and the energy at which the 
prompt emission spectrum peaks, £p ea k (Amati et al. 2002), we 
conclude that this relation, and the prompt emission in general, 
does not provide a reliable measure of the total energy. As a 
corollary, we rule out the narrow jet opening angles used by 
Lamb, Donaghy & Graziani (2004), Oj ~ 0. 1° and thus confirm 
that the true GRB rate is significantly lower than the rate of type 
Ib/c SNe (Berger et al. 2003a). 

Finally, the overall agreement between the energies derived 
here and those inferred from modeling of the relativistic phase 
of the afterglow indicates that the central engine in GRBs 
980703 and 970508 did not produce a significant amount of en- 
ergy in mildly relativistic ejecta (T/3 J> 2) at late time, t ~ ?nr- 
However, a comparison to the beaming-corrected 7-ray ener- 
gies (Bloom, Frail & Kulkarni 2003), £ 7 » 1.1 x 10 51 erg 
(GRB 980703) and £ 7 - 10 51 erg (GRB 970508) reveals that 
the efficiency of the blastwave in producing 7-rays, e 7 , varies 
considerably: - 20% for GRB 980703, but only - 3% for 
GRB 970508. The wide dispersion in e 7 strengthens the con- 
clusion that £ 7 is not a reliable tracer of the total energy (Berger 
etal. 2003b). 

The low value of e 7 for GRB 970508 may indicate an in- 
jection of energy from mildly relativistic ejecta at early time. 
Both the optical and X-ray light curves of this burst exhibited 
a sharp increase in flux approximately 1 day after the burst, by 
a factor of about 4 and ;> 2, respectively (Piro et al. 1998; 
Sokolov et al. 1998). The flux in these bands depends on en- 
ergy as F v oc E (p+3> / 4 and oc E (p+2) I A , respectively (Sari, Piran 
& Narayan 1998). Thus, if we interpret the flux increase as due 
to injection of energy from ejecta with V ~ 5- 10 (Panaitescu, 
Meszaros & Rees 1998) we find an energy increase of about 
a factor of three. The analysis performed here provides an es- 
timate of the total energy following the injection and thus e 7 
appears to be low. The actual value of e 7 is thus ~ 10%. 

Although GRBs 980703 and 970508 are currently the only 
bursts with sufficient radio data to warrant the full Sedov- 
Taylor analysis, flattening of radio light curves at late time have 
been noted in several other cases, most notably GRBs 980329, 
991208, 000301C, 000418 and 000926 (Frail et al. 2004). In- 
terpreting the flattening as a transition to non-relativistic expan- 
sion and using the expression for the flux at 8.5 GHz at the time 

of the transition, F v {tm) ~ 50[(1 +z) /2]~ l l 2 t e -ieg\n i J A Esid^ 
/iJy (Livio & Waxman 2000), we find the rough results 

n Q /4 E 5l « 6 (980329), w 4 (991208), « 25 (00030 1C), « 6 
(000418), and w 22 (000926). Thus, for typical densities, 
~ 1 - 10 cm -3 (Panaitescu & Kumar 2002; Yost et al. 2003), 
the inferred kinetic energies are again of the order of 10 51 - 10 52 
erg. 

This leads to the following conclusions. First, the energy 
scale of cosmological bursts is about 5 x 10 51 erg, at least three 
orders of magnitude higher than the kinetic energies in fast 
ejecta determined for local type Ib/c SNe from radio observa- 
tions (Berger, Kulkarni & Chevalier 2002; Berger et al. 2003a), 
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and an order of magnitude higher relative to the nearby (d w 40 
Mpc) GRB 980425 associated with SN 1998bw (Kulkarni et al. 
1998; Li & Chevalier 1999; Waxman 2004b) and GRB 031203 
(z = 0.105; Prochaska & et al. 2004; Soderberg & et al. 2004). 
Second, as already noted in the case of GRB 030329 (Berger et 
al. 2003b), there is a wide dispersion in the fraction of energy 
in ultra-relativistic ejecta, such that the 7-rays are a poor proxy 
for the total energy produced by the engine. 

Thus, radio calorimetry is uniquely suited for addressing 
the relation between various cosmic explosions. So far, such 
studies reveal a common energy scale in relativistic ejecta of 
about 5 foe (foe = 10 51 erg) for cosmological GRBs (Berger et 
al. 2003b), about 0.1 foe for the low redshift bursts (980425, 



031203), and £ 10" J foe in fast ejecta for type Ib/c SNe. The 
open question now is whether we are beginning to trace a con- 
tinuum in the energetics of cosmic explosions, or whether the 
various classes truly represent distinct physical mechanisms 
with different energy scales. Fortunately, the best example to 
date of an object possibly bridging the various populations, 
GRB 030329, still shines brightly in the radio a year after the 
burst. 

We thank Eli Waxman, Sarah Yost and Re'em Sari for valu- 
able discussions, and the referee, Roger Chevalier, for useful 
comments. We acknowledge NSF and NASA grants for sup- 
port. 
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Table 1 

Physical Parameters of GRBs 980703 and 970508 



Parameter 


GRB 980703 


GRB 970508 


r(10 17 cm) 


1.05-2.5 


3.7-5.9 


B(G) 


0.02-0.7 


0.04-0.25 


7 


8-270 


65-165 


n (cm -3 ) 


8-3.5 X 10 3 


0.4-10 


e e 


0.01-0.45 


0.07-0.5 




5 x Kr 6 -0.4 


1 x nr 3 -o.45 


Mej,i so (10" 4 M S ) 


1-40 


2-18 


£ S T (10 50 erg) 


9-56 


15-38 


£x(?dec) (10 51 erg) 


4 


30 



NOTE.— Physical parameters of GRBs 980703 and 970508 derived 
from the non-relativistic evolution of their blastwaves. The range of al- 
lowed radii, and hence physical parameters, is determined by the condition 
(E e +Eb) < Est/2. The last entry in the table, £x(fdec). is the total kinetic 
energy at the deceleration time, « dcc 90 s, including synchrotron radiative 
losses (§5). 
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io 1 10 2 10 3 

Time Since the Burst (days) 

Fig. 1.— Radio light curves of the afterglow of GRB 980703 at 1.4, 4.9 and 8.5 GHz. Only data at t > r NR = 40 d (black circles) are 
used in the fit. The data exhibit a clear flattening relative to the relativistic evolution of the afterglow (thin gray line; Frail et al. 2003) 
in agreement with the expected change from F v oc t~ p (jet) to F v oc t ( 2l ~ l5 P)/ m (ISM) or F v oc t ( - 5 ~ 7p ^ 6 (Wind) in the non-relativistic 
regime. The best-fit light curves for the ISM (black) and Wind (gray) models are indistinguishable. The models include a contribution 
from the host galaxy of 40, 50 and 65 /iJy at 8.5, 4.9 and 1.4 GHz, respectively. 
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0.5 1 1.5 2 2.5 3 

Radius (cm) x 10 17 

Fig. 2.— Energies associated with the afterglow of GRB 980703 in the non-relativistic Sedov-Taylor phase as a function of the 
(unconstrained) blastwave radius. The thin curve is the sum of the energy in relativistic electron (E e oc r~ 6 ) and in the magnetic fields 
(Eb oc r 11 ). Also plotted are the Sedov-Taylor energy (£st oc r~ 2 ) and the thermal component, Est/2. The shading corresponds to 
an uncertainty of 30% in the value of the synchrotron frequency u at t = r NR . The value of E ST /2 provides an additional constraint, 
E e + E B < Est:/ 2, which limits the range of allowed radii in the solution (boxed region). 
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Fig. 3.— Physical parameters of the Sedov-Taylor blastwave for GRB 980703 at f NR = 40 d for the range of radii that obey the 
constraint E e + Es < £st/2 (Figure 2): (a) The ratio of energy in the relativistic electrons to that in the magnetic fields, (b) the 
magnetic field strength, (c) the density of the circumburst medium, (d) the Sedov-Taylor energy, (e) the velocity of the blastwave, 
and (f) the isotropic-equivalent mass of the ejecta produced by the central engine and responsible for the afterglow emission. The 
light shaded region in (c) marks the range of densities inferred from the relativistic evolution of the fireball, n w 20-35 cm" 3 (Frail 
et al. 2003). With the additional constraint that the density derived here conform to this value, we derive the values of e e /e B , B, Est, 
(3, and M e j marked by arrows. 
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Fig. 4.— Radio light curves of the afterglow of GRB 970508 at 1.4, 4.9 and 8.5 GHz. Only data at t > r NR = 100 d (black circles) 
are used in the fit. The best-fit light curves for the ISM model are shown (black); the Wind model can be ruled out since it requires 
p<2. 
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Radius (cm) x10 17 

Fig. 5. — Energies associated with the afterglow of GRB 970508 in the non-relativistic Sedov-Taylor phase as a function of the 
(unconstrained) blastwave radius. The thin curve is the sum of the energy in relativistic electron (E e oc r~ 6 ) and in the magnetic fields 
(Eg oc r 11 ). Also plotted are the Sedov-Taylor energy (Est oc r" 2 ) and the thermal component, Est/2. The shading corresponds to 
an uncertainty of 30% in the value of the synchrotron frequency uo at t = r NR . The value of Est/2 provides an additional constraint, 
E e + E B < Est/2, which limits the range of allowed radii in the solution (boxed region). Finally, the arrow marks the most likely 
solution using the value of the cooling frequency as estimated from a combination of the radio and optical data (§4). This additional 
parameter breaks the radius degeneracy, indicating r w 4.2 x 10 17 cm and £ S x ~ 3 x 10 51 erg 
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Fig. 6.— Physical parameters of the Sedov-Taylor blastwave for GRB 970508 at ?nr = 100 d for the range of radii that obey the 
constraint E e + E B < E ST /2 (Figure 5): (a) The ratio of energy in the relativistic electrons to that in the magnetic fields, (b) the 
magnetic field strength, (c) the density of the circumburst medium, (d) the Sedov-Taylor energy, (e) the velocity of the blastwave, 
and (f) the isotropic-equivalent mass of the ejecta produced by the central engine and responsible for the afterglow emission. The 
arrows mark the most likely values using an estimate of the cooling frequency from a combination of the radio and optical data (§4 
and Figure 5). 



